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We show that the density matrix of a spin-Z system can be described entirely in terms of the 
measurement statistics of projective spin measurements along a minimum of 4/ + 1 different spin 
directions. It is thus possible to represent the complete quantum statistics of any iV-level system 
within the spherically symmetric three dimensional space defined by the spin vector. An explicit 
method for reconstructing the density matrix of a spin-1 system from the measurement statistics of 
five non-orthogonal spin directions is presented and the generalization to spin-i systems is discussed. 
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I. INTRODUCTION 



As rapid progress is being made in the experimental generation of quantum states, it becomes necessary to develop 
efficient methods of characterizing the actual mixed state output of each new realization. In particular, various types 
of optical spin-1 systems have recently been generated using parametric downconversion 0, 0, 0, [J IE @ • It is therefore 
interesting to consider the measurements necessary to properly identify the quantum states of such spin systems. 

In the most general case, these states can be characterized by reconstructing the complete density matrix from a 
sufficiently large set of measurements, a procedure commonly referred to as quantum tomography fj\. 111. Igl. Il0|. For 
two- level systems (qubits) , quantum tomography is usually realized by measuring the three orthogonal components of 
the Bloch vector represented by the Pauli operators. In spin-1/2 systems, the physical meaning of these components is 
generally clear. In particular, they represent the components of the three dimensional Stokes vector in the commonly 
studied case of single photon polarization 0, ■ In spin-? systems with higher total spin, the connection between the 
much larger number of density matrix elements and the physical properties of the system is less clear. For abstract 
N-level systems (qudits) , an expansion of the density matrix into the generators most closely related to the individual 
density matrix elements has been proposed §J. However, the physical properties corresponding to these operators are 
quite different from the spin components observed e.g. in Stern-Gerlach or n-photon polarization measurements. 

In particular, the recently generated n-photon polarization states are usually characterized by photon detection 
measurements in a pair of orthogonal polarization directions QLftldQ- This corresponds to a projective measurement 
of one component of the three dimensional Stokes vector, which is formally equivalent to the three dimensional spin 
vector. The direction of the Stokes vector component is determined by the pair of orthogonal polarization directions 
detected in the measurement and can be varied by using standard birefringent elements such as quarter wave plates 
and half wave plates. The experimental characterization of such optical spin-/ systems thus corresponds to the 
■ measurement of spin components Li along a set of well defined measurement directions i. It is therefore desirable to 
formulate quantum tomography in terms of the measurement statistics obtained in this kind of measurements. 

In the following, we show how the measurement statistics obtained in measurements of spin components Li relate 
to the elements of the density matrix. Based on these results, a systematic approach to the quantum tomography 
of spin-Z systems is developed. We propose a decomposition of the density matrix into components that reflect the 
spherical symmetry of the spin system and correspond directly to well-defined contributions in the experimentally 
observable spin statistics. It is shown that measurements along a minimum of Al + 1 spin directions are necessary 
to reconstruct the complete density matrix. An explicit method for reconstructing a spin-1 density matrix from the 
measurement probabilities along five non-orthogonal spin directions is derived and the extension of this method to 
general spin-Z systems is discussed. Since this method can be applied equally well to small (few level) and large 
(many level) quantum systems, it also provides a useful basis for the discussion of decoherence and the transition 
from quantum to classical physics. 
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II. MEASUREMENT STATISTICS OF A SPIN-Z SYSTEM 

Each projective von-Neumann measurement of a spin component Li results in one of the 21 + 1 eigenvalues m of 
the quantized spin along the direction corresponding to i. By repeating the measurement a large number of times, 
it is possible to determine the probability distribution Pi(m) of the measurement outcomes m. The information 
represented by this probability distribution can also be expressed in terms of averages of different powers of L i: 

duT) = j2 m " M™)- (!) 

m 

The probability distribution over the 21 + 1 possible outcomes is then uniquely defined by the 21 averages obtained 
for n = 1 to n = 21, that is, the 21 expectation values form a set of linearly independent parameters describing the 
complete measurement statistics. 

Using this representation of the measurement information obtained along one spin direction, it is now possible to 
derive the relations between the measurement statistics along different spin directions by expressing the spin statistics 
of an arbitrary spin direction in terms of the three orthogonal spin components L x , L y , and L z . If the spin direction 
is given in terms of the horizontal and azimuthal angles <f> and 6, the measurement statistics of L((f>, 9) then read 

^ (L(<t>, °)) n ) = ( (sin(0) cosO)!^ + sin(0) sm.(<f>)L y + cos((9)£ 2 ) . (2) 

Each nth order expectation value ((Li) n ) can therefore be expanded into expectation values of nth order products 
of the three orthogonal spin components. Specifically, the linear spin expectation values (n = 1) are defined by the 
three components of the average spin, (L x ), (L y ), and (L z ), while the quadratic spin expectation values are given 
by expectation values such as (L x L y + L y L x ), describing correlations and fluctuations of the spin, and so on. The 
measurement statistics of any measurement direction can therefore be described by a hierarchy of expectation values 
ranging from n = 1 to n = 21. 

Although there are an infinite number of measurement directions, only a finite number of independent parameters 
is necessary to describe all nth order expectation values. In general, these independent parameters can be obtained 
by multiplying out equation J3J) and identifying the contributions associated with different functions of the angles <p 
and 6. Alternatively, it is possible to identify these contributions with elements of the density matrix. As will be 
described in more detail below, this corresponds to an expansion of the density matrix into basis operators A„^ that 
can be expressed as nth order polynomial functions of the spin component operators Lj. Since the density matrix 
uniquely determines the spin statistics, the number of independent nth order expectation values should be equal to 
the number of independent basis operators A„^ in the expansion of the density matrix. The additional density matrix 
elements associated with an increase of spin from spin-(Z — 1/2) to spin-/ can then be identified with the additional 
(n = 2l)th order expectation values required to describe the spin statistics given by equation (j5J). 

The number of parameters defining the density matrix of any iV-level system is N 2 — 1. A spin-Z system is an 
N = 21 + 1 level system, and the number of parameters defining the density matrix is 41(1 + 1). A spin-(/ — 1/2) 
system is an N — 21 level system, and the number of parameters defining the density matrix is 4l 2 — 1. The number 
of additional parameters needed to describe a spin-Z system rather than a spin- (7 — 1/2) system is therefore equal 
to 4/ + 1. Since these additional parameters characterize the (n = 2Z)th order of the measurement statistics given 
by equation ©, the number of independent nth order expectation values should be equal to 2n + 1. Note that this 
result is indeed consistent with the number of symmetrically ordered nth order products of the three spin components 
required to describe (£"), once the constant value of the total spin length defined by L 2 + L 2 + L 2 = 1(1 + 1) has 
been taken into account. Explicit examples for components of the nth order spin expectation values will be given in 
sections IIIII and 

Table [I] illustrates the distribution of density matrix parameters. Any spin-Z system is characterized by 2n + 1 
expectation values of nth order, with n running from n = 1 to n = 21. That is, the density matrix can always be 
represented by three linear spin averages, five second order spin averages, seven third order spin averages, and so 
on [HI. Experimentally, each measurement along a given spin direction determines one nth order average for each 
value of n. Since complete quantum tomography requires the determination of 41 + 1 independent contributions to the 
expectation values of order n = 21, it is therefore necessary to measure at least 41 + 1 different spin directions in order to 
obtain the necessary measurement information |l2j . Note that this condition is a result of allowing only measurements 
of spin components Li. If general von-Neumann measurements were possible, 21 + 2 measurements would be sufficient 
for quantum tomography. In the case of spin component measurements, about half of the information obtained is 
redundant, since it reproduces the results for lower order expectation values that can be obtained from fewer spin 
directions. For example, quantum tomography of a spin-2 system requires measurements along nine spin directions, 
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Components of the 




System size 


density matrix 






(L-) (T?) (T, 3 ) (f^) ■■• 


total 


spin-1/2 (2-level) 


3 _ _ _ ... 


3 


spin- 1 (3-1gv61) 


3 5 


8 


spin-3/2 (4-level) 


3 5 7 


15 


spin- 2 (5-level) 


3 5 7 9 


24 



TABLE I: Reconstruction of the density matrix from n-th order expectation values of the spin obtained from projective 
measurements of various spin directions i. The total number of matrix elements necessary for tomography is iV 2 — 1, where 
N — 21 + 1 is the dimensionality of the corresponding Hilbert space. 

providing nine averages to determine the three linear spin expectation values, nine averages to determine the five 
quadratic spin expectation values, nine averages to determine the seven third order spin expectation values and nine 
averages to determine the nine fourth order expectation values. Thus it is only the need to determine the complete 
(n = 2Z)th order spin statistics that makes it necessary to measure a total of Al + 1 spin directions Lj. 

Note that the precise choice of measurement directions is not very critical, since the only requirement for obtaining 
the complete spin statistics is that 2n + 1 of the 4/ + 1 averages obtained from the nth order statistics are linearly 
independent. In some special cases, the information obtained from Al + 1 different measurement directions may not 
be sufficient to reconstruct the complete density matrix because the choice of measurement directions coincides with 
a symmetry in the (n = 2Z)th order expectation values. For example, the 2nd order information obtained for L z is 
already obtained from measurements of L x and L y , since L x + L 2 + L 2 = 1(1 + 1). However, even a small tilt of 
one of the measurement axes can fix this problem by providing additional information. Therefore, almost any choice 
of Al + 1 measurement directions will be sufficient for quantum to mog raphy. To minimize errors, it may be useful 
to keep the angles between different directions as large as possible However, the choice of Al + 1 measurement 
directions used for a complete reconstruction of the density matrix can generally be quite arbitrary, and each specific 
tomography protocol merely represents one example out of an infinity of equally valid possibilities. 

In order to obtain an explicit description of the density matrix in terms of nth order spin statistics, it is necessary 
to identify the contributions of different order n in the density matrix. This can be achieved by expanding the density 
matrix using an appropriate operator basis. In general, there are infinitely many expansions of the density matrix 
into orthonormal basis operators. For reasons of mathematical simplicity, the most common choice is that given by 
the generators of the SU(iV) algebra [H|l3- 

However, these generators do not represent the spherical symmetry of 
the spin-Z system. We therefore propose an alternative expansion of the spin-Z density matrix p\ that is based on the 
different orders n of the spin statistics, given by basis operators X n ,i, where n represents the lowest order of the spin 
statistics to which X n ^ contributes, and i is the index of the component within this order, running from 1 to 2n + 1 
for each value of n. The conditions for orthogonality and for the normalization of this basis then read 

Tr{A nii } = 

Tr{A n; iA m j} = *2 5 n rn 5ij . (3) 

Using this complete operator basis, it is possible to expand the spin-Z density matrix pi in terms of expectation values 
of the spin statistics as 

21 /2n+l \ 

A = 2zTl i +2E E<Vi>An,i • (4) 

n=l \ i=l / 

It is thus possible to formulate the density matrix entirely in terms of the measurement statistics of the three di- 
mensional spin vector, reflecting the analysis of the measurement statistics described by equation given in table 
[I] Specifically, the (./V 2 — l)-dimensional basis is divided into groups of nth order products of the spin operators with 
2n+ 1 operators each, describing the separate nth order contributions to the spin statistics. In the following, we show 
how such a basis of three linear and five quadratic spin operators can be defined for the specific case of I = 1. 
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III. EXPANSION OF THE SPIN-1 DENSITY MATRIX 



The case of I = 1 provides the most simple example of non- linear contributions to the spin statistics. Moreover, 
optical spin-1 systems have already been realized by two photon polarization states 0, 0, OJ or by the orbital 
angular momentum of single photons The following procedure for quantum tomography of a spin-1 system may 

therefore be particularly useful in the characterization of such experimental results. 

For 1 = 1, the linear spin components already fulfill the conditions defined by equations The remaining five 
basis operators can then be constructed using 2nd-order operator products. One such set of quadratic operators that 
fulfills the relations defined by 10 is 15] 
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In this basis, the generalized eight dimensional Bloch vector can be separated into a three dimensional linear part 
and a five dimensional quadratic part, given by 
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In the L 2 -basis, the matrix elements of these operators read 
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(7) 



The expectation values of these eight operators characterize the density matrix in terms of the linear and the quadratic 
measurement statistics given by equation . The linear expectation values along any spin direction are given by the 
three expectation values of Ai j, 



{L{cj>, 6)) = sin(0) cos(0)(A M ) + sin(0) sm(^)(\ h2 ) + cob(0)<Ai, 3 }, 
and the quadratic expectation values are given by the five expectation values of X2A, 



(8) 



(L{^ef 



+ - sin 2 ((9) cos(2<?!>)(A 2il ) + sin 2 (0) sin(0) cos(</»)(A 2 , 2 > + sin(6>) cos(6») sin(0)(A 2 , 3 ) 
+ sin(0) cos(0) cos(0)(A 2 , 4 ) + -J=(l - | sin 2 (0))(A 2 , 5 ). 



(9) 
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Using this relation, it is possible to determine the correct expectation values of all five 2nd order basis operators 
A2.i from the quadratic expectation values of five independent measurement directions. Together with the linear 
expectation values, the results of these five measurements then define the complete density matrix, 

1 1 3 1 5 

Pi=i = o 1 + o £>i,i>Ai,i + 2 I>2,i>A2,< • (10) 



from (L(<t>,e)) from (L(<f>,e) 2 ) 

An explicit procedure for quantum tomography can now be formulated by choosing a set of five measurement directions. 
The expectation values of the basis operators X nt i can then be expressed in terms of the measurement probabilities 
along the five measurement directions. 



L\ — L x , Li — L y , L% — —p^ ( L x + L y 

l 4 = ^(4 + 4), L 5 = j=(L z + i 



IV. QUANTUM TOMOGRAPHY OF THE SPIN-1 SYSTEM BASED ON THE MEASUREMENT 

STATISTICS OF FIVE SPIN DIRECTIONS 

For 1 = 1, each measurement along a given spin direction Li has three possible outcomes, m = ±1 and m = 0. In 
the experimentally relevant case of two-photon polarization [ij, |2j, y, |4j , these measurement outcomes correspond to 
the detection of two horizontally polarized photons (m = +1), two vertically polarized photons (m = +1), and one 
photon each in horizontal and in vertical polarization (m = 0), where the component of the Stokes vector is selected 
by appropriate rotations of the polarization using standard linear optics elements. The measurement statistics given 
by the probabilities Pi(m) can then be identified with normalized coincident count rates at the detectors. 

As discussed above, five measurement settings are necessary to perform quantum tomography. A particularly simple 
choice of the five measurement directions for the spin-1 system is given by 

-4 ( l x + Lj] , 

(11) 

In the case of two-photon polarization, the first three spin directions can be identified with linear polarization rotated 
by zero for L±, by 7r/4 for L 2 , and by 7r/8 for L3. The remaining two directions then represent elliptical polarizations, 
with the main axes along ±7r/4 for £4, and the main axes along angles of zero and ir/2 for i 5 . Using this set of 
measurement settings, it is now possible to explicitly identify the expectation values (A n> i) that define the density 
matrix with the corresponding measurement probabilities pi(m). 

In general, the measurement probabilities of a spin-1 system along a given direction i are related to the nth order 
expectation values of the corresponding spin component Li by 

(Li) = Pi(+1)- Pi (-1) 

(4 2 ) = l-ft(0). (12) 

The quadratic terms of the spin statistics are thus entirely determined by the measurement probabilities Pi(0) for 
a measurement value of zero spin along the measurement direction i. The five expectation values A2,i defining the 
quadratic components of the density matrix can therefore be obtained from the five measurement probabilities Pi{0) 
along the spin directions i = 1 to i = 5 defined by i|llfl . The relations between the measurement statistics and the 
expectation values of the corresponding basis operators then reads 

-( Pl (0)-p 2 (0)) 
Pi(0)+p 2 (0)-2p 3 (0) 
Pi(0)-2p4(0) + l 
p 2 (0)-2p 5 (0) + l 
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(13) 



The measurement probabilities for m — thus determine five of the eight coefficients in the expansion of the density 
matrix given by equation Q). The remaining three coefficients can be obtained from the linear expectation values 
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given by the differences between Pi(+l) and 1), e.g. 
(Ai,i) = (L x ) = P1 (+1)_ P1 (-1) 
(Ai, 2 ) = (4) = P2(+1)-Pa(-1) 

(A 1)3 ) = (i») - -V2((P3(+1)-J» 3 (-1))-(P4(+1)-P4(-1))-(P5(+1)-P5(-1))). (14) 

Since the results of five measurement directions are used to determine only three parameters, there are two relations 
between the measurement results that should be approximately fulfilled if the measurement error is low. These 
relations can be written as 

( Pl ( + 1) « -L ((P3(+1) - pa(-l)) - (p 4 (+l) - P4 (-l)) + (P6(+l) - P5(-l))) 

(P2(+1)-P 2 (-1)) « -^((P3(+1)-P3(-1)) + (P4(+1)-P4(-1))-(P6(+1)-P6(-1))). (15) 

Effectively, the measurement results for the spin directions three to five can be used to predict the results for spin 
directions one and two. Relations (|15|l thus illustrate the application of quantum tomography to the prediction 
of further measurements on the same system. The experimental differences between this prediction and the actual 
outcome of the measurements may therefore provide a realistic estimate of the errors limiting the reliability of quantum 
tomography in practical applications. 

Since the eight parameters (\ n ,i) completely define the density matrix, they can be used as an alternative rep- 
resentation of the quantum state, just like the three dimensional Bloch vector for two level systems. Any other 
representation of the density matrix can then be obtained from equation if the representations of each element of 

the basis is known. For example, the density matrix elements in the L z -basis given by equations J7J) can be used to 

express any density matrix defined by the expectation values (\ n ,i) as 

pt=i = 

5 + 3<*i.8> + 575<*a.«> 5^ «A M > - i(A lia ) - i<A3, 3 > + <Aa,4» § «A 2 ,i) - i(A 2 , 2 » 

-^((A lil )+i(A 1 , 2 )+i<A 2 , 3 ) + (A 2 ,4)) 1-^(^,5) ^ 5 ((A 1 ,i)-*<Ai,2)+ l <A 2 ,3>-(A 2 ,4» . 

§«A 2 ,i)+*<A 2 , 2 » _^((A 1 , 1 ) + l (A 1 , 2 )-i(A 2 , 3 )-(A 2 , 4 )) 5-5^.3) + 575<W) 

(16) 

The quantum coherences between different eigenstates of L z can thus be expressed in terms of averages involving 
the other two spin components, L x and L y . In fact, this relationship between quantum coherence in the L z -basis and 
the orthogonal spin components L x and L y can be used to systematically construct density matrix decompositions 
for higher spins, as we will show in the following. 



V. QUANTUM STATISTICS FOR I > 1 

In order to generalize the construction of a convenient set of basis operators X n ^ to arbitrarily large spins, it is 
useful to organize the basis operators according to their density matrix elements in the L 2 -basis. Such an organization 
is already indicated by the example for spin-1 given in equation In this example, each group of density matrix 

elements | m)(m' | with the same order of coherence \m! — m\ in L z depends on a well defined subset of the operators 
\ ny i- This organization of basis operators can be generalized to arbitrary spins I by formulating the operator basis 
for each spin value in such a way that the matrix elements of each operator are non-zero for only one value of 
\m! — m\. The relationship between the hierarchy of spin expectation values (\ n ,i) and the density matrix elements 
in the L z -basis is then as shown in table [n] Specifically, the 21 + 1 diagonal matrix elements (\m! — m\ = 0) are 
determined by the set of 21 + 1 operator expectation values (A n ,(2n+i)) that can be obtained from the measurement 
statistics of L z according to equation Q). Likewise, the 41 off-diagonal elements with \m! — m\ — 1 can be determined 
by the set of 41 operator expectation values, (X n ,2n) and (A^n-i), and so on. Note that the 41 — 2 off-diagonal 
elements with \m — m\ = 2 do not include any first order spin statistics, since the linear spin components only have 
matrix elements up to \m' — m\ = 1. Likewise, the 41 — 4 off-diagonal elements with \m! — m\ = 3 do not include 
any linear or quadratic spin statistics. In general, an expectation value of at least |m' — ro|th order is necessary to 
describe the effect of a corresponding off-diagonal element in the spin statistics. 
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TABLE II: Relation between the density matrix elements (m' | p \ m) in the L z -basis and the corresponding nth-order basis 
operators representing the spin statistics. 



A convenient way to construct the two operators A n l and A nj 2 with non-zero density matrix elements of maximal 
coherence \m — m!\ = n is to apply the non-hermitian spin operator L x + iL y . This operator only has non-zero matrix 
elements with mf — m = +1. Specifically, 

[L x + it y ^ | m) = - m)(l + m + 1) | m + 1). (17) 

Consequently, (L x + iL y ) n has only matrix elements with m! — m = n. It is therefore possible to generate the nth 
order basis operators A Mi i and X n ^ with matrix elements of \m' — m\ = n from the normalized hermitian components 
of {L x +iL y ) n , 



An 1 — 



\L X + iLy) n + (L x iLy) n 
\Jtt{{L x + iL y ) n (L x - iL v ) n } 
-i[[L x +iL y ) n -(L x -iL y ) n ^ 

y/Tr{(L x + iL y ) n (L x -iL y ) n } 



(18) 



Starting from these definitions of basis operators, the complete set of basis operators may be constructed, e.g. by 
multiplying the operators A n !/ 2 with different powers of L z and/or L x + L y to obtain higher order contributions with 
the same coherence \m' — m\ in the density matrix. The precise factors can be determined using the requirements 
for orthogonality and normalization given by equation (jSJ). It is then possible to construct a complete orthonormal 
operator basis for any spin-/ system. 

By establishing the relation between spin statistics and coherence in the density matrix, equation 118fl also illustrates 
the physical meaning of quantum coherence in L z . In particular, it is worth noting that the greater the difference 
\m' — m\ between the L z - eigenvalues of the states that are in a coherent superposition, the higher the order of the 
spin expectation values in L x and L y that is needed to identify this coherence in the measurement statistics. The spin 
correlation hierarchy presented in tabel[n]may thus provide a key to understanding the non-classical effects associated 
with quantum superpositions in arbitrarily large physical systems. 



VI. NON-CLASSICAL CORRELATIONS AND DECOHERENCE 



The highest possible value for \m! — m\ in the density matrix is obtained for quantum coherence between the extremal 
.£,z-eigenstates | m — +1) and | m = —I). According to tabel|nl the two matrix elements describing this coherence 
correspond to the two (2/)th order spin expectation values. These expectation values can be constructed explicitly 
using equation l|18l) and include only products of L x and L y . In order to characterize a coherent superposition of 
the L 2 -eigenstates | m = +1) and | m = —I), it is therefore necessary to evaluate the (2Z)th order spin statistics 
in the xy-plane. For all orders lower than 21, the spin statistics obtained in the measurements of such a coherent 
superposition are identical to those of an incoherent mixture of | m = +/) and \m = —I). 
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This observation has significant implications for the identification of strong non-classical effects in large quantum 
systems. At sufficiently high values of I, the superposition of | m = +1) and | m = —I) is a cat-state like superposition 
of two macroscopically distinguishable states. It is therefore interesting to know that the effects of this superposition 
appear only in the highest order expectation value of the spin statistics. While the lower order expectation values are 
very easy to measure since only very few measurement directions are required and the measurement errors tend to 
average out, the highest order expectation values can only be determined from sufficiently precise measurement results 
of at least 41 + 1 measurement directions. Effectively, the highest order expectation values represent a measurement 
resolution at the quantum level, providing the information necessary to resolve the precise eigenvalues of the spin 
components 0|. This means that the (2Z)th order spin statistics is very sensitive to errors of ±1 in the spin 
measurements. In other words, the smallest measurement errors are sufficient to make the effects of the cat-state like 
superposition between | m = +1) and | m = —I) disappear. We can therefore conclude that the actual non-classical 
properties of a superposition of macroscopically distinguishable states can only be observed in the microscopic details 
of the measurement statistics. It is therefore not surprising that decoherence quickly wipes out such tiny details. 

For a more precise evaluation of decoherence and measurement precision, it is useful to consider the case of isotropic 
decoherence caused by spin diffusion due to random rotations. The time evolution of the density matrix caused by 
this kind of decoherence can be described by 

j t Pi = -T Yl {I^P+^-UplX (19) 

i=x,y,z 

Using the well known commutation relations of the spin operators, it is possible to calculate the relaxation dynamics 
of the nth-order expectation values of the spin. For the non-Hcrmitian operators (L x + iL y ) n the result reads 



^LitT , .? \n\ _ n(n+l) 
dt 



-((L x + iL y ) n ) = — L_^r((£ lt + iL y r). (20) 



Since the relaxation of the spin is isotropic, all nth order contributions to the expansion of the density matrix should 
relax at the same rate. The effect of isotropic decoherence therefore reduces each nth order parameter (X n ,i) of the 
density matrix expansion Q by a decoherence factor of exp[— Ttn(n + l)/2], and the time evolution of the density 
matrix can be written as 

1 21 T n(n+ 1)1 ( 2n+1 \ 

m = 2J+1 1 + ^ exp ~ Tt 2 ^ (*»,i>t=o Ki • (21) 

n=l <- -J \ i=l / 

The expansion of the density matrix pi into an operator basis X n ^ based on the different orders of the spin statistics 
therefore greatly simplifies the description of any isotropic errors in the preparation and manipulation of spin states. 

Since the decoherence effects described by equation (|21ll arise from spin diffusion, it is also possible to identify Tt 
with an increasing uncertainty in the spin direction, Tt = S9 2 /2. The result of equation (|21|l can then be used to 
estimate the errors caused by a misalignment of the measurement direction. Specifically, an alignment error of 59 
will reduce the expectation values observed for the nth order spin statistics by a factor of exp[— 56 2 n(n + l)/4]. To 
obtain at least exp[— 0.25] = 78% of the original expectation value at orders n 3> 1 of the spin statistics, the errors of 
the spin alignment have to be smaller than 59 — 1/n. The precision in the alignment of the spin direction necessary 
to obtain the nth order statistics is thus proportional to 1/n and the requirement for observing evidence of a cat-like 
superpositions in spin-/ systems is an angular resolution of 59 < 1/(21). 



VII. ENTANGLEMENT STATISTICS AND GENERAL SPIN NETWORKS 



The formalism developed above can also be applied to entangled spin-/ systems. In this case, the density matrix 
of the total system is obtained by evaluating the correlations between measurements of the local spin components. 
Specifically, the joint quantum state of a spin-Z^ system A and a spin-Z^ system B can be determined by simultaneously 
measuring spin components Li(A) in A and spin components Li(B) in B, obtaining the joint probabilities pij(uiA, tub) 
of each measurement outcome. The correlated spin statistics can then be expressed in terms of the expectation values 

((Li(A)) riA (Li(A)) nB ) = m T ™ n B B Pij{mA,m B )- (22) 

m.A,m B 

It is then possible to analyze the spin statistics according to the local order tia and ub, where the total number of 
independent components required to characterize each order is given by the product (2n,4 + l)(2ns + 1). Note that 
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in this case, ua = and tlb = have to be included in order to describe the local spin statistics of each system. 
Consequently, the lowest order expectation values are given by (ua — = 0) and (tia = 0,tib = 1), with three 

independent components each. The second order expectation values (tia + ns = 2) are given by five components for 
(tia = 2, tib = 0), nine components for (tia = 1, tib — 1), and five components for (tia = 0, ns = 2). The highest order 
contribution to the correlated spin statistics is then given by (ha = 21a, ns — 21b), with a total of (41a + 1)(4/b + 1) 
independent components. The number of measurement settings required to perform complete quantum tomography 
for entangled spin systems is therefore equal to (41a + l)(4Zs + 1). In the experimentally realized case of I a = Ib = 1 
00, 5, 6], this would require 25 different measurement settings with nine possible outcomes each, for a total of 225 
measurement probabilities. 

An explicit description of the density matrix in terms of the correlated (nA,ns)t\\. order spin statistics can be 
obtained using products of the basis operators for each individual system. The expansion of the density matrix then 
reads 

pAB = Wa + mi B + D i ® i + 2WTT) g — ® l) x — ® l ) 

+ 2(2/ + i) S 8 8 + 4 12 \J2CK A ,i^k B ,j)K A ,i®KBd • (23) 

The expectation values defining the density matrix can now be expressed in terms of the joint measurement proba- 
bilities •pij( m A,' m 'B) by writing the (71^4, ri^ )th order expectation values of the correlated spins in equation J22J) as a 
function of the expectation values in equation Q23JI • It is then possible to fully characterize any NxM entanglement in 
terms of the correlated spin statistics. 

The extension of this formalism to multi-partite spin networks is also straightforward, since the density matrix 
can be expanded into products of the local basis operators for any number of systems. The expectation values of 
these products can then be determined from the correlated measurement statistics of spin measurements performed 
simultaneously on all systems. 



VIII. CONCLUSIONS 



In conclusion, we have shown how the density matrix of spin-Z systems can be reconstructed from the measurement 
statistics of projective spin measurements along a set of at least 41 + 1 different spin directions. The components of 
the density matrix can then be identified with different contributions to the statistics of the three dimensional spin 
vector. It is therefore possible to interpret the discrete quantum statistics of arbitrarily large spin systems within the 
same three dimensional space defined by the Bloch vector of a two level system. 

The explicit procedure for the quantum tomography of spin-1 systems provides an example of the general method 
that can be applied directly to experimentally generated two photon polarization states such as the ones reported in 
0, 0, El- EI ■ It may thus serve as the foundation of a more detailed characterization of decoherence and noise effects 
in these newly available entanglement sources. 

Besides its practical usefulness for the experimental characterization of general spin-Z systems, the expansion of the 
density matrix into elements of the spin statistics also provides a more intuitive understanding of quantum statistics 
in large systems. The analysis presented above may therefore also help to clarify the conditions for the emergence of 
quantum effects in physical systems of arbitrary size. 

Acknowledgements 

H.F.H. would like to thank A.G. White and K. Tsujino for some very motivating discussions. 



[1] T. Tsegaye, J. Soderholm, M. Atature, A. Trifonov, G. Bjork, A. V. Sergienko, B. E. A. Saleh, and M. C. Teich, Phys. 
Rev. Lett. 85, 5013 (2000). 

[2] A. V. Burlakov, M. V. Chekhova, O. A. Karabutova, and S. P. Kulik, Phys. Rev. A 64, 041803(R) (2001). 

[3] A. Lamas-Linares, J.C. Howell, and D. Bouwmeester, Nature 412, 887 (2001). 

[4] J. C. Howell, A. Lamas-Linares, and D. Bouwmeester, Phys. Rev. Lett. 88, 030401 (2002). 



10 



[5] A. Mair, A. Vaziri, G. Weihs, and A. Zeilinger, Nature 412, 313 (2001). 
[6] A. Vaziri, G. Weihs, and A. Zeilinger, Phys. Rev. Lett. 89, 240401 (2002). 

[7] A.G. White, D.F.V. James, P.H. Eberhard, and P.G. Kwiat. Phys. Rev. Lett. 83, 3103 (1999). 
[8] D.F.V. James, P.G. Kwiat, W. J. Munro, and A.G. White, Phys. Rev. A 64, 052312 (2001). 
[9] R.T. Thew, K. Nemoto, A.G. White, and W. J. Munro, Phys. Rev. A 66, 012303 (2002). 
[10] G. Klose, G. Smith, and P.S. Jessen, Phys. Rev. Lett. 86, 4721 (2001). 

[11] It may be worth noting that, in terms of group theory, the 2n + 1 expectation values of nth order form an irreducible 
representation of the spherical symmetry defined by the arbitrary direction of the spin component Li in equation 10 . 

[12] It has recently come to our attention that a particular choice of 4/ + 1 measurement directions was proposed as early as 
1968 by R. G. Newton and B.-L. Young, Ann. Phys. (N.Y.) 49, 393 (1968). The results in this paper show that 41 + 1 
measurement directions are also sufficient for quantum tomography. However, Newton and Young suggested that a better 
choice of directions could reduce the number of measurement directions needed. Our results show that this is in fact 
impossible. 

[13] F.T. Hioe and J. H. Eberly, Phys. Rev. left. 47, 838 (1981). 

[14] G. Mahler and V. A. Weberrufi, Quantum Networks (Springer, Berlin 1998), p. 44. 
[15] H.F. Hofmann, Phys. Rev. A. 68, 034307 (2003). 

[16] In fact, the truncation of the measurement statistics at {(Li) 21 } itself is a consequence of the fact that the discreteness of 
the quantized eigenvalues makes a higher resolution unnecessary. 



